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ABSTRACT 

We reconsider correlations among the spectral peak energy {Ep), 1-second 
pe ak luminosity (Lp ) and isotropic energy (i^iso ), using the database constructed 



by lYonetoku et al.l ( 120101 ) which consists of 109 Gamma-Ray Bursts (GRBs) 



whose redshifts are known and Ep, Lp and Ei^o are well determined. We divide 
the events into two groups by their data quality. One (gold data set) consists 
of GRBs with peak energies determined by the Band model with four free 
parameters. On the other hand, GRBs in the other group (bronze data set) have 
relatively poor energy spectra so that their peak energies were determined by the 
Band model with fixed spectral index (i.e. three free parameters) or by the Cut- 
off power law (GPL) model with three free parameters. Using only the gold data 
set we found the intrinsic dispersion in logLp (= cxint) is 0.13 and 0.22 for Ep-T^- 
Lp correlation (Tl = E-^so/Lp) and Ep-Lp correlation, respectively. We also find 
that GRBs in the bronze data set have systematically larger Ep than expected by 
the correlations constructed with the gold data set. This means that the intrinsic 
dispersion of correlations among Ep, Lp, and i^iso of GRBs depends on the 
quality of data set. At present, using Ep-T^-Lp correlation with gold data set, 
we would be able to determine the luminosity distance with ~ 16% error, which 
might be useful to determine the nature of the dark energy at high redshift 2 > 3. 
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Subject headings: (cosmology:) distance scale — (stars:) gamma-ray burst: gen- 
eral 



Introduction 



Discoveries of empirical correlations of gamma-ray bursts (GRBs) raised many re- 
searches on early universe using GRBs. One of the most well studied correlations is the 
one between the sp ectral peak energy (E^) and isotropic equivalent energy ( E\^n ) called 



Ep-E\f,n cor r elation (lAmati et al.ll2002l : ISakamoto et al.ll2004j : iLamb et al.ll2004j : lAmati et al. 



20061 . 12OO9I ). lYonetoku et al.l (120041 ) found a similar but tighter correlation between Ep and 



1-second peak luminosity called the Ep-Lp correlation. These correlations are tight but they 
have large dispersions such as aint = 0.33 in logLp and ajnt = 0.37 in logE 'i sn which can not 
be explained as statistical errors of Ep , Eiso and Lp (lYonetoku et al.ll2010l ). iGhirlanda et al. 



( I2OO4] ) found that E p tight ly correlates with the coUimation-corrected gamma-ray energy 
(E^). iFirmani et al.l (120061 ) proposed that adding the high signal time scale (T0.45) to the 
Ep-Lp relation reduces the dispersion of the correlation. This correlation is defined by us- 
ing only prompt emission properties like Ep-Eiso , Ep-Lp correlations so that it seems to 
be promising tools to constrai n the cosmologica l parameters. However, th is correlation is 
not confirmed by lat er studies (IRossi et al.ll2008l : ICoUazzi fc Schaeferll2008l ). More recently, 
Tsutsui et al.l (120091 ) found that the luminosity time (Tl = Eiso/Lp) also improves both the 



Ep-Eis 



and Ep-Lp correlations. 



These correlations were us ed to investigate the st ar formation history ( lYonetoku et al. 



2OO4J ). the reionaization ep och (iMurakami et al 
tory of the early universe (iTakahashi et al.l 12003 



20051). and the cosmological expansion his 



tory 01 the early universe (llaKahashi et ai.MzUUdi: lUguri 
2006 : ISchaefeJ2007l : lKodama et al.ll2008l : lLiang et al.lboosi : 

mm- 



Qguri fc Takahashi 



20061: 



Cardone et al.l 



Ghirlanda et al. 



200 



Tsutsui et al 



However, in spite of high correlation coef ficients, there have been many cautions to 



use these erapirical correlations for cosmology (iNakar fc PiranI |2004J : iBand &: Preecd 12005 



Butler et al.ll2007l : IShahmoradi fc Nemiroflal2009l ). To establish these correlations in GRBs 
prompt emissions as tools to determine cosmological parameters, we must investigate the 
origins of systematic errors and the way to remove them. We note that there are many 
factors to cause systematic errors besides intrinsic dispersions of their prompt emissions. 
For example the sensitivity of the detectors , the evolution effects of GRBs, the confusion 
with other sources, the lack of unknown parameters like the jet opening angle djet, etc. All 
of these effects might arise the additional systematic errors over the intrinsic dispersions of 
GRBs. 
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Possible se lection effects on t 
trast i ng results (Butler et al.lboOTl: 



lese correlations are studied by many authors with con- 



2009: Amati et al. 



12009 



Ghirlanda et al.ll2008l : iNava et al.ll2008l : IShahmoradi fc Nemiroff 
Yonetoku et al. 201ol ). However previous st udies did not conside r 



the difference of spectral models to determine Ep well. As shown in iKaneko et al.l ( 120061 ). 
it often happens that high energy power-law index /3 for the Band model with four free 
parameters can not be determined by the data so that the cutoff power-law (GPL) model 
with three free parameters is used to fit the data. GPL model might be good if the peak 
energy is close to the high energy end of the detector band width. One can not use Band 
model but GPL model if the event is so dim that t he nu mber of high energy photons is very 
small. Importantly, simulations in iKaneko et al.l (120061 ) showed that, if the signal-to-noise 
ratio is relatively low , a true spectrum with the shape of the Band model can be fitted by 

100 keV. Therefore 



GPL model with E"""' 



which is larger than the true value of E^^^ up to 



GPL model might overestimate Ep. While , if we fit a true GPL spectrum by the Band 
model, the estimated value of is almost equ al to the true value since the la rge value of 
— /3 looks like an exponential function. In reality IShahmoradi fc NemirofB ( l2009l ) found that 



Ep estimated using the GPL model by IKaneko et al.l (120061) are systematically harder than 



Ep estimated using the Band model by lYonetoku et al.l ( 12004 ). Although the systematic 
difference between the peak energies fitted by the Band model and the ones fitted by the 
GPL model are reported, how this difference affect the spectral-brightness correlations of 
GRBs has been hardly studied so that we shall study this problem in this paper. 

The purpose of this letter is to investigate the effect of uncertainty in using different 
spectral models which determine on the Ep-E^f^q , Ep-Lp and Ep-T^-Lp correlations. 



using our database developed in lYonetoku et al.l (120101 ). We examine this model bias by 
dividing the samples into two data sets as gold and bronze according to the quality of 
spectral observation. In this paper, we assume, if signal-to-noise ratio is high enough, all of 
the spectrum of GRBs are well expressed by the Band function. 

The structure of this letter is as follows. First we describe our database of 109 GRBs 
with known redshift and well-determined spectral parameters, 1-second peak luminosity, and 
isotropic energy in section |2j We construct the Ep-T^-Lp , Ep-Lp 



and E„-E^^r, correlations 



with only gold data set in section |31 Finally we will give summary in section |H 



2. Data Description 



In lYonetoku et al.l (120101). we constructed a database selecting 1 09 GRBs from GGN 
Gircular Archive (iBarthelmyl Il997l ) and GRBlog (iQuimby et al.ll2003l ). In this section we 
briefly describe our database. 
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Let us begin with Ep. In many cases, the prompt gamma-ray spectrum is well fitted 
wi th the spectral mo del of the exponentially-connected broken power-law function suggested 
by lBand et al.l (119931 ) . This Band function has four parameters, the low-energy photon index 
a, the high-energy photon index P, the spectral break energy Eq and the normalization A. 
The peak energy (Ep), at which the flux is maximum in the uF^ spectrum, can be calculated 
as Ep = {2 + a)Eo. 

However, for some GRBs, the photon index (mostly /3) cannot be deter mined due to the 



limited energy range of the detector and/or the lack of the number photons ( jPendleton et al. 



19971 ). When the observation of high-energy range is not enough, the spectrum is sometimes 



fitted with the Cut-off power law (CPL) function. This function has three parameters, the 
low-energy photon index a, the spectral break energy Eq and the normalization A. In this 
case the peak energy can be derived as Ep = (2 + a)EQ. Note that even if, for a given 
GRB spectrum, the reduced chi square value of this model is smaller than that of the Band 
function, it is difficult to say whether this model reflects the intrinsic property of the GRB 
or it is just due to the poor statistics in the high-energy range. The reported values of 
Ep for GRBs which were poorly observed in the high-energy range are based on either the 
Band function or CPL function, depending on the observation team so that there exists the 
ambiguity in the definition of Ep from the biginning. 



Let us move on to Eiso, Lp and Tl. In lYonetoku et al.l (120101 ). we calculated the bolomet- 
ric energy and the peak luminosity in the energy range 1-10,000 keV in the rest frame of each 
GRB by extending the observed spectrum. Here, it should be noted that the integration was 
performed assuming the Band function even for GRBs whose spectra were not fitted by the 
Band function and the photon indices were not reported. In these cases we assumed the typi- 
cal values a = —1 and f3 = —2.25 to calculate the bolometric fiuence (Shr, \) and the bolomet- 
ric peak flux Fp boi- These values are suggested by BATSE observations (jPreece et al.ll2000[ ) 
and also supported by Fermi observations of GRB 080916C, 081024B, 090323 and 090428 
up to possibly 100 GeV energy range. Then the bolometric isotropic energy (-Eiso) and the 
1-second peak luminosity (Lp) can be simply calculated as i^iso = 47rc?|5'boi/(l + z) (erg), 
and Lp = 47r(i|Fp boi (erg s~^). Here, cIl is the luminosity distance calculated for the fiat uni- 
verse with the cosmological parameters of {flm, ^a) = (0.3, 0.7) and the Hubble parameter 
of Hq = 70 km s~^Mpc~^. Further we define the luminosity time as the third parameter of 
GRB prompt emission as Tl = Eiso/Lp. The error of the luminosity time is estimated by 
using error propagation equation. We can neglect the crossterm b etween Lp and . Eiso b ecause 
of the independence of the Ep-E-^^o and Ep-Lp relation shown in (ITsutsui et al.ll2009l ). 

Thus, for GRB whose observed photon number is small, there are two possible systematic 
effects. One comes from the fact that the peak energy Ep is determined by fitting the 
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spectrum with either the Band function or CPL function. As iKaneko et al.l (|2006[ ) pointed 
out that the CPL function tends to overestimate Ep compared to the Band function. This 
would induce a systematic error in the correlations related to Ep. On the other hand, 
although Lp and £^iso are determined in a single straightforward way, the photon indices are 
set to the typical values if the number of detected photons is small. This would also cause 
a systematic error. 

To estimate these systematic errors, we call a certain GRB belongs to the gold data 
set if its spectrum is well observed so that it is fitted by the Band function quite well and 
all four parameters are accurately determined. Other GRBs for which the fixed a and /or 



are allocated to bronze data set. Here 3 of 109 GRBs in the database of lYonetoku et al. 



( 120101 ). do not have i^iso and then Ti so they are included in neither the gold nor the bronze 
data sets. As a result, the number of the gold and bronze data sets are 41 and 65 GRBs, 
respectively. 

In the following sections, we construct the Ep-Eiso , Ep-Lp and Ep-TL-Lp relations for 
the gold data set. It is expected that the correlations obtained from the gold data set would 
suffer from relatively small systematic errors so that we could study real intrinsic dispersions 
of the correlations. 



3. correlations 

Here we derive the Ep-T^-Lp , Ep-Lp , and Ep-Ei^o correlations from the gold data set. 

First we assume the correlation among Ep, and Lp to be of the form, logLp = 
A + i?log (i?p/440 keV) + Clog (Ti/4.70 s), where we take the denominator of the second 
and third terms as the average value of Ep and Tl in gold data set to minimize the statistical 
errors for these correlations. For this correlation, we found seven outliers (980425, 980613, 
000131, 090328, 091003, 091020, 091127) shown by blue color in top left of Fig.l. The seven 
outliers deviate from the best-fit relation at more than 3-cr dispersion level. We will give 
some arguments on this point in section HI Using 34 gold data set of GRBs, we calculate the 
best fit function shown by the solid black line in top left of Fig. 1 with red points of the gold 
data set and 3 a errors of the Ep-T^-Lp relation by the yellow color region. The functional 
form of the best fit function is given by 

/ p \ l-70±0.07 f rp \ -0.40±0.06 

L„ = (52.64 ±0.03) X ( ^\ ( tit^ I ■ (1) 

^ ^ ' V440keVy V4-70secy ^ ' 

Here, we include not only errors in Lp but also errors in Ep and so that the chi-square func- 
tion is defined as x^{A B, C) = S(logL°^^ - A - 5 log (^p/440 keV) - Clog {Tl/4.70 s))V(f^mcas+ 
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crfnt) where the first term of weighting factor is ericas = (1 + '^C)'^\ogLp + (BaiogEp) + 
(CfTiogTi)^- The factor 2C in front of cfogL comes from the fact that the definition of 
includes Lp. The reduced chi-square is unity wit h the intrinsic dispe rsion o"int = 0.13. This 



correlation is consistent with our previous study iTsutsui et al.l (|2009[ ). 



Similarly, we can obtain the best-fit function and errors of the Ep-Lp (top right of Fig.l) 
and Ep-E[so (bottom of Fig.l) correlations for the same 34 gold GRBs, 



1.76±0.10 



= (52.63 ± 0.05) X (^^) , (2) 



with the intrinsic dispersion (Tint = 0.22, and 



1.68±0.13 



= (53.31 ± 0.06) x(^^ J , (3) 

with the intrinsic dispersion (Tint = 0.31. The best fit values of the Ep-Lp and Ep-Eiso 
correlations are consistent w i th pre vious studies but the intrinsic dispersions are tighter 



than those in lYonetoku et al.l ( 120101 ) in which both gold and bronze data are used. 



The values of (Tint suggest that Ep-T^-Lp correlation is tightest among three correlations. 
In Fig. 1, we show the gold data set (red points) with the best-fit function (solid line) and 
3-sigma dispersion region (dotted lines). We see that the Ep-T^-Lp correlation is much 
tighter than the Ep-Lp and Ep-E\so correlations by eye also. The blue points and green 
points indicate seven outliers of the Ep-T^-Lp correlation and bronze data set, respectively. 
The bronze data set are systematically harder and/or dimmer than the gold data set. This 
difference causes a large dispersion in addition to the intrinsic dispersion of the correlations 
if we include the bronze data in the analysis. 



4. Summary &z Discussion 



In this paper, using database constructed by lYonetoku et al.l (l2010l ). we examine the 
model bias, that is. Band or GPL, on Ep-Ti-Lp , Ep-Lp and Ep-E\so correlations. We found 
that GRBs with the peak energies fitted by the GPL model are distributed in systematically 
harder and/or dimmer side of the Ep-T^-Lp , Ep-Lp and Ep-Eiso correlations than the ones 
by the Band function. There might be two interpretations about this result. The first is that 
these correlations have much larger intrinsic dispersion than that of observed one. If we had 
the more sensitive de tector and could observe dimmer GRBs, the disp ersion of the relations 
would become larger (iButler et al.ll2007l : IShahmoradi &: Nemirofal2009l ). Another is that the 
use of the GPL model to estimate the peak energies causes this systematic difference. As 
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Fig. 1. — The Ep-T^-Lp relation (top left), the Ep-Lp correlation (top right) and the Ep-E\so 
(bottom) correlation with all data set. The solid line and dotted lines indicate the best fit 
function and 3-0" dispersion region in Eq ([T]), Eq ([2]) and Eq ([3]). The bronze data set seems 
to be harder and/or dimmer than gold data set. The CPL model or the Band function with 
fixed power-law index cause this systematic difference. 
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simulated by iKaneko et al.l (|2006[ ). the Band function spectrum is well fitted by the CPL 
model if the detector does not have enough sensitivity to observe the high-energy photons. 
However the peak energies fitted by the CPL mo dels are always high er than that of the 
simulated Band function spectrum (see table. 3 in iKaneko et al.ll2006l ). Thus, it seems to 
be natural to conclude that the latter is more acceptable. In short, using only the peak 
energies determined by the Band function, we would get tighter correlations. If we could 
have much more GRBs by which we can uniformly analyze the data with the Band function, 
GRBs would be more powerful tool to constrain cosmological parameters. 

We found seven outliers in our gold data set. We classify these outliers in two classes as 



(dimmer and/or harder) 980425, 980613, 090328, 091003 
(brighter and/or softer) 000131, 091020, 091127 



Although we do not know how and why these outliers are different from ordinary GRBs 
except for the distribution in the Ep-T^-Lp space, the effect of eliminating these GRBs is 
obvious, that is, the correlation becomes tighter. To find the characteristics which distinguish 
these outliers from ordinary GRBs is urgent. We here point out possible origins of these 
outliers. Let us assume that if we observe the jet nearly on axis , Ep-T^-Lp correlation 
would be very tight. However if we observe the jet with a certain viewing angle, we might 
have some dispersions on the observed Ep-T^-Lp correlation. In other words, we might 
not avoid some dispersion in the Ep-T^-Lp correlation from viewing angle, especially when 
the observer locates near the edge of the jet of GRBs. If we will know how to distinguish 
these outliers from the ordinary gold GRBs, the Ep-T^-Lp correlation might be much tighter 
and very useful in determining the nature of dark energy in redshift larger than ~3. We 
should note t hat even in the P eriod-Luminosity relation of Cepheid variable there are ~ 
10% outliers ( iRiess et al.l 120091 ) so it is not surprising that there are ~ 20% outliers in the 



Ep-Ti-Lp relation. 



Butler et al.l (120071 ). using the Bayesian approach to estimate Ep, indicated that dim 
events close to the detector sensitivity would make large scatter on the Ep-Eis^ and Ep-Lp 
relations and that there is a significant threshold effect. Thus, they conclude that the Ep-Ei^o 
correlation have larger intrinsic dispersion thari observed if we do not suffer from a threshold 
effect. Recently, IShahmoradi fc Nemirofla (120091 ) argue that using hardness ratio instead of Ep 
they also find that £'p'"*-Fluence correlation become more wider if we will be able to determine 
Ep of dimmer events. However, there is a possible bias by using different method to estimate 
Ep. Even the difference of Ep between the Band and CPL models causes the systematic 
errors so that using the other method to estimate Ep might cause the additional systematic 
errors. The smaller the intrinsic dispersion of the relation is, the more the correlation suffers 
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from these syste matic effects. This might be why the Ep-Tn,4;^-Lp relation is not confirmed 
by later studies (IFirmani et al.l 120061 : iRossi et al .1120081 : ICoUazzi fc Schaeferll2008[ ). 



Kaneko et al. (2006) suggested the Ep value of CPL function becomes systematically 
higher than the one of the Band function. If the 65 bronze data previously analyzed by the 
CPL function are reconsidered by the Band function with the fixed (3 as an average value 
of —2.25 (Preece et al. 2000), they might show the distribution around the best fit line of 
each correlation estimated with 41 gold data set. They have a good potential to become a 
"silver" data set. To do so, we need help from each instrument team, and this is a future 
work. 

Finally, we note that there would be many reasons which cause systematic errors on 
the correlation in addition to intrinsic property of GRBs. These systematic errors must be 
carefully estimated and removed from the correlation analysis one by one. If we will finish 
it, GRBs become more powerful and unique standard candles to investigate the nature of 
the dark energy at high redshift larger than ~ 3. 
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